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Abstract 

We consider the stability of synchronized chaos in coupled map lattices and in 
coupled ordinary differential equations. Applying the theory of Hermitian and 
positive semidefinite matrices we prove two results that give simple bounds on 
coupling strengths which ensure the stability of synchronized chaos. Previous 
results in this area involving particular coupling schemes (e.g. global coupling 
and nearest neighbor diffusive coupling) are included as special cases of the 
present work. 

PACS Numbers: 05.40.+j, 05.45.+b 
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I. INTRODUCTION 



Synchronization of coupled chaotic systems has found applications in a variety of 
fields including communications optics f|||, neural networks |||| and geophysics 

||. An essential prerequisite for these applications is to know the bounds on the coupling 
strengths so that the stability of the synchronous state is ensured. Previous attempts [|10| -|29| 
aimed at obtaining such conditions have typically looked either at systems of very small 
size or at very specific coupling schemes (diffusive coupling, global all to all coupling etc. 
with a single coupling strength). More recently, Pecora and Carroll introduced the 
notion of a master stability function for general coupling topologies, but this function can 
only be accessed in a numerical fashion. The contribution of this Letter is a methodology 
that can lead to analytical bounds on the individual coupling strengths under the mild 
assumption that the coupling be symmetric. We demonstrate the method in the form of 
two theorems, one applicable to coupled map lattices and the other to coupled ordinary 
differential equations (ODEs). 



II. STABILITY RESULTS 



The coupled map lattices we consider are in the form 

1 L 



x>+l) = f(x>)) + f(x*(n))-f(x«(n)) , (1) 

L 3=1 



and we use 



± i (t)=f(x i (t))+j:b ij *(t), (2) 

3=1 

for coupled ODEs. Here x* is an M-dimensional state vector describing the ith map/ODE 
and f is an M-dimensional map/vector field. We assume that the coupling is symmetric, 
and bij = bji, and that, in the absence of coupling, the individual M-dimensional 
system is chaotic with the largest Lyapunov exponent h max > 0. With the additional 
constraint that H,=i hj = it is easy to see that the synchronized state, x 1 = x 2 = • • • = 
x L = x, is a solution to our models. We obtain stability conditions by requiring that the 
"transverse" Lyapunov exponents (which will be defined later) are all negative. We note 
that the stability obtained here is in the sense of weak stability (or Milnor stability), which 
ensures the transverse stability of typical orbits, and is a necessary condition for asymptotic 
stability where all orbits are transversely stable. Our stability results for the synchronized 
chaotic solution are expressed in the following two theorems. 

Theorem 1. For Eq. (|l|) the synchronized chaotic state is stable if for all i, j, j ^ i, 

[1 - exp(-h max )) < aij < [1 + exp(-h max )]. (3) 

Theorem 2. For Eq. (0) the state of synchronized chaos is stable if 

kj > h max /L, V i, j, j ^ i. (4) 
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We now sketch the proof of Theorem 1. Linearizing Eq. (HD around the synchronized 
chaotic state x we get : 



z i (n + 1) = J(x(n)) 



(5) 



where z* denotes the ith map's deviations from x and J is the M x M Jacobian matrix. In 
terms of the M x L state matrix S(n) = (z 1 (n) z 2 (n) . . . z L (n)) Eq. ([|) can be written as 
the following matrix equation: 

S(n+1) = J(x(n))S(n)C T , (6) 
where C T is the transpose of the LxL coupling matrix C containing the coupling coefficients: 

[C]« = 1 - -J^Oy, l=l,2,...,L, 

[C]y = aij/L, i^j (7) 

Since the coupling coefficients are assumed to be real and symmetric, C is Hermitian. 
It can be diagonalized and all its eigenvalues are real C = EAE -1 . Here E and A are 
the eigenvector and eigenvalue matrices of C, respectively. Let e be one of the eigenvectors 
of E and A its associated eigenvalue. Acting Eq. (0) on e we get: 



Let u(n) = S(n)e. Then 



S(n + l)e = AJ(x(n))S(n)e. (8) 



u(n + l) = AJ(x(n))u(n). (9) 



We now compute the Lyapunov exponents for the above reduced system. We note that 
A = 1 is always an eigenvalue of C and its corresponding eigenvector is (1 1 ... 1) T . In this 
case, the above equation is just the linearization of the individual map which was assumed 
to be chaotic. Therefore, the eigenvector (11 ... 1) T of C with eigenvalue 1 corresponds to 
the synchronized chaotic state. The Lyapunov exponents in this case are nothing but the 
Lyapunov exponents for the individual system. Hence they are given by hi = h max , h 2 , ... , 
Hm- These describe the dynamics within the synchronization manifold defined by x* = x Vi. 

Next we consider the remaining eigenvalues and eigenvectors. Since C is a symmetric 
matrix, the remaining eigenvectors span a (L — l)-dimensional subspace orthogonal to the 
eigenvector (1 1 . . . 1) T . Consequently, this subspace is orthogonal to the synchronization 
manifold. For each A ^ 1 we calculate the Lyapunov exponents for Eq. ([5]). Since A is a 
real number (C being a symmetric matrix), the Lyapunov exponents are easily calculated. 
Denoting them by /xi(A), /U 2 (A), /ijvf(A), we have 

Hi(X) = hi + \n\X\, z = l,2,...,M. (10) 



We refer to these Lyapunov exponents as transversal Lyapunov exponents [O since they 
characterize the behavior of infinitesimal vectors transversal to the synchronization manifold. 
These determine the stability of the synchronized chaotic state. For stability, we require the 



3 



transversal Lyapunov exponents for each A 7^ 1 to be negative. This is equivalent to the 
statement 



/Wr( A) = h max + In I A| < 0. (11) 

In other words, we require |A| < exp(— h max ) for each A 7^ 1. From this equation we see that 
for the synchronized chaotic state to be stable, A = 1 should be the eigenvalue of C with 
the largest magnitude. Ordering the eigenvalues of C as Ai = 1 > A 2 > A 3 > • • • > A^, the 
stability conditions can be rewritten as 

A 2 < exp(-/wr), (12) 
Al > — exp( h max ). (13) 

Our goal is to obtain bounds on such that the above two inequalities are simultane- 
ously satisfied. Consider a Hermitian matrix K defined as follows: 

[K] U = 1- {L ~ L 1)R , z = l,2,...,L, 

[K]« = f. ( 14 ) 

where R is a constant which will be characterized later. Consider the matrix P = K — C. 
We see that the diagonal elements are given by 

P« = ^E««- - f = l,2,...,L. (15) 

These are positive if > R V i, j, j 7^ i. Next consider the absolute value of the off-diagonal 
elements of P: 

|P*;i = lf-^l, Vui/i. (16) 



If ciij > R, it can be seen that 



|P«|>EI P «I> ^ = l,2,...,L. (17) 



This implies that P is positive semidefinite . 

We now introduce the concept of positive semidefinite ordering Since Hermitian 

matrices are generalizations of real numbers and positive semidefinite matrices are gener- 
alizations of nonnegative real numbers, one can introduce an ordering among Hermitian 
matrices as follows |JT]: Let A, B be L x L Hermitian matrices. We write A y B if the 
matrix A — B is positive semidefinite. Further, if A >z B, a\ > 02 > ■ ■ ■ > oil are the 
ordered eigenvalues of A and Pi > fa > • • • > Pl are the ordered eigenvalues of B, then 

(Xi>Pi, i = l,2,...,L. (18) 

Earlier, we had already shown that P = K — C is positive semidefinite if > i? V i, j j 7^ 
i. Since K and C are Hermitian, we have K y C. Now the largest eigenvalues of both K 
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and C are equal to 1. The second largest eigenvalue of K can be easily calculated and is 
equal to 1 — R. Therefore K y C implies that 

\ 2 <{1-R). (19) 

Comparing with the inequality given in Eq. (|12|), we see that this constraint is obeyed if 
(1 - R) < exp(-h max ). That is, 

R>l-exp(-h max ). (20) 

But A 2 < 1 — R only if > R V i, j j ^ i. Putting the two inequalities together, we get 
the first stability condition: 

Uij > [1 - exp(-h max )], V z, j, j ^ i. (21) 

Next we need to satisfy the second stability constraint given in Eq. (0). Following a 
procedure similar to the one above, we get the second stability condition: 



a 



it 



< [1 + exp(-/i ma2 )], V j ^ i. (22) 



Combining the stability conditions given in Eqs. fl2T|) and (|22|) we get our final result as 
follows. The synchronized chaotic state of Eq. ([I]) is stable (in the Milnor sense) if the 
coupling coefficients (which are assumed to be symmetric) obey the stability condition 
given in Eq. (H). Note that as h max — > 0, the range of allowed coupling strengths increases 
reaching a maximum range of (0,2). As h max — > oo, the range decreases to zero. This 
confirms the intuitive expectation that it should be harder to stabilize a synchronized state 
which is more chaotic. The above result generalizes earlier results [0,1151 obtained assuming 



that all coupling coefficients are identical. Further, in some of the earlier papers ||11|| , it 
was assumed that the constant (dimensionless) coupling strength is less than 1 and hence 
the upper bound given here (which is greater than 1) was not explicitly observed in those 
previous studies. 

We now numerically verify the above result by studying a system of 100 Henon maps. 
The / [cf. Eq. ([[])] for an individual Henon map is given by: 

fi(x 1 ,x 2 ) = 1 + x 2 - ax\] f2(xi,x 2 ) = bx 2 . (23) 

For the values a = 1.4 and b = 0.3, the maximum Lyapunov exponent h max is found to be 
0.43. We now couple 100 Henon maps using the scheme given in Eq. (|l]) where we randomly 
generate the coupling strengths a^-'s. From Eq. (|3|), the synchronized chaotic state is stable 
if 0.35 < ciij < 1.65. We have numerically verified this result. 

If we have nearest neighbor coupling we can obviously obtain better bounds than given 
in Eq. (H) since we know more information about the coupling coefficients. In this case, 
we obtain the following bounds (further details on this and other coupling schemes can be 
found in our forthcoming paper [[3^]): 



1 - exTp(-h max ) l + exp(-h max ) 

— — — t—— < an < — — — if Li is odd, 

1-cos(2tt/L) 1 1 + cos(2?r/L) ' 

1 - exp(-h ma x) l + exp(-h max ) . 

; -— < clu < , if L is even. 

1 - cos(2vr/L) 3 2 
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Note that as L becomes larger, the above range shrinks to zero. We find [|32| a conservative 
estimate for the critical value L c which makes the range zero to be: 

mt (cos- 1 [exp(-h max )/ (1 + exp(-/i mtw .))] ) ' ^ 
This generalizes the result found in Ref . [JTTJ] . The above result implies that for a sufficiently 



large L (> L c ) nearest neighbor coupled systems can not have a stable synchronized chaotic 
state. 

We now prove Theorem 2. The initial treatment is similar to the one used by Pecora and 
Carroll in arriving at the master stability function |3(|. The structure of coupling that we 



have assumed includes the commonly used diffusive coupling, nearest neighbor coupling, all- 
to-all coupling, star coupling etc. Our proof breaks down if only one or few of the components 
of x l are coupled. In this case, Pecora and Carroll PDA have shown numerically that more 



complicated stability conditions arise. Note that, unlike the coupled map case, we only have 
a lower bound in the stability condition. This difference arises from the fact that it is only 
for maps the stability condition is in terms of the absolute value of the eigenvalue which 
leads to both lower and upper bounds. This is not so for coupled differential equations. 

The proof of this theorem is along the same lines as our proof for coupled maps. Lin- 
earizing around the synchronized state we get 

z* = J(x)z l + 5> jZ J , (25) 

where z l denotes deviations from x and J is the M x M Jacobian matrix. We now introduce 
the M x L state matrix S = (z 1 z 2 ... z L ). Then the linearized equation Eq. ( 25|) can be 
written as the following matrix equation: 

S = J(x)S + SC T , (26) 

where C T is the transpose of the LxL coupling matrix C containing the coupling coefficients: 

[C]u = — ^ bij, i = 1, 2, . . . , L, 

[C]ij = In-, 1 i ± j (27) 

Since the coupling is symmetric, C is Hermitian. Following the same procedure as before, 
we obtain 

u = [J(x) + AI] e, (28) 

where I is the MxM identity matrix, e is an eigenvector of C and A its associated eigenvalue. 
We now compute the Lyapunov exponents for the above reduced system. We note that 
A = is always an eigenvalue of C and its corresponding eigenvector is (1 1 ... 1) T . This 
corresponds to the synchronized chaotic state. The Lyapunov exponents in this case are 
nothing but the Lyapunov exponents for the individual system. Hence they are given by 
hi = h maxj h 2 , ... , hjrf. 

Next we consider the remaining eigenvalues and eigenvectors. Since C is a symmetric 
matrix, the remaining eigenvectors span a (L — l)-dimensional subspace orthogonal to the 
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synchronization manifold. For each A 7^ we calculate the Lyapunov exponents for Eq. (f28|). 
Since XI commutes with J(x), the Lyapunov exponents are easily calculated. Denoting them 
by /ii(A), // 2 (A), /U M (A), we have 

fi i (X) = h i + X, z = l,2,...,M. (29) 

For stability, we require these transversal Lyapunov exponents for each A 7^ to be negative. 
This is equivalent to the statement 

Hrnax(X) = h max + A < 0. (30) 

In other words, we require A < —h max for each A ^ 0. Ordering the eigenvalues of C as 
Ai = > A 2 > A 3 > • • • > Xl, the stability conditions can be rewritten as 

A 2 < -h max ■ 

(31) 

As before, we wish to obtain bounds on 6^- such that the above inequality is satisfied. 
Consider a Hermitian matrix K' defined as follows: 

[K.% = -{L-l)R\ * = 1,2,...,L, 

[K% = R>, i^j, (32) 

where R' is a constant which will be characterized later. Consider the matrix P' = K' — C. 
We see that the diagonal elements are given by 

P« = £fci i=l,2,...,L. (33) 

These are positive if > R' V i, j, j 7^ i. Next consider the absolute value of the off-diagonal 
elements of P': 

|Py = |#-M, (34) 
If bij > R', it can be shown that 

|Py>£l p y> i = l,2,...,L (35) 



and P' u > 0. This implies that P' is positive semidefinite [ pT| . 

We have shown that P' = K' — C is positive semidefinite if b^ > R' V i, j j 7^ i. Since 
K' and C are Hermitian, we have K' y C. Now the largest eigenvalues of both K' and C 
are equal to 0. The second largest eigenvalue of K' can be easily calculated and is equal to 
—LR'. Therefore K' y C implies that 

A 2 < -LR'. (36) 

Comparing with the inequality given in Eq. (0), we see that this constraint is obeyed if 
-LR' < -h max . That is, 

R > h max /L. (37) 

But A 2 < —LR' only if by > R! V i,j j 7^ i. Putting the two inequalities together, we get 
the required stability condition given in Eq. (01). 
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III. CONCLUSIONS 



To conclude, we have derived very simple bounds on the coupling coefficients which 
ensure the stability of the synchronized chaotic state of L symmetrically coupled systems. 
We also gave specific bounds for the nearest neighbor coupled map system. These results 
allow for non equal coupling coefficients and generalize earlier results found in the literature 



ilj,13|l which were obtained assuming that all coupling coefficients are constant. It is very 



easy to apply our criteria to the system being studied and they encompasses a wide class 
of coupling schemes including most of the popularly used ones in the literature. Further, 
we expect the introduction of non-equal coupling to lead to interesting new phenomena 
in coupled systems. Our stability results would enable a systematic exploration of such 
systems. 

Our results were made possible by a sequence of operations. We summarize them below 
since we feel that this approach is applicable to the stability analysis of a wide variety of 
coupled systems and not merely the specific problem considered in this letter. First, we 
converted the linearized system to a matrix equation. This was further simplified by looking 
at the evolution of eigenmodes of the coupling matrix C. By realizing that only the largest 
Lyapunov exponent matters for our analysis, the stability conditions for the synchronized 
chaotic state were recast as bounds on certain eigenvalues of the coupling matrix. Then we 
bound the Hermitian matrix C in the coupled map system and the coupled oscillator system 
by carefully constructed constant Hermitian matrices K and K' respectively. This was done 
in such a manner that that K — C and K' — C are positive semidefinite when the coupling 
coefficients satisfy certain inequalities. Using a powerful result from matrix analysis, this 
automatically implies that the eigenvalues of C are bounded above by the eigenvalues of K 
(or K'). By comparing the different bounds that we had derived, we then finally arrived at 
Eqs. (0) and (^) which ensure the stability of the synchronized chaotic state in the coupled 
map and coupled oscillator systems respectively. 
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